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Abstract 

We study the dynamics of particles coupled to gravity in (2 + 1) dimensions. Using the ADM formalism, 
we derive the general Hamiltonian for an iV-body system and analyze the dynamics of a two-particle system. 
Nonlinear terms are found up to second order in k in the general case, and to every order in the quasi-static 
limit. 

1 Introduction 

The study of 2+1 dimensional gravity [I] [2] 3 4 5, 6, 7 8 ( 9] has been around for several decades. Its attraction 
was rooted in its mathematical simplicity, which afforded some insight into the dynamics and behaviour of 
general relativistic gravity. Indeed, in this framework, the Einstein tensor can be expressed in terms of the 
curvature tensor, such that vacuum must be locally flat. This framework became even more popular after the 
discovery of the BTZ black hole [101 111] and is now a standard tool practitioners of quantum gravity employ in 
understanding their subject [121 1131 PH] . 

We are concerned here with the problem of iV-body dynamics, which is a long-standing one in relativity 
due to its notorious difficulty. In lower dimensional settings, this problem becomes much simpler. For example, 
the general form of the solution has been obtained for lineal gravity [TS] for arbitrary N, after which a variety 
of exact solutions for N = 2 were obtained in various contexts that include both charge and cosmological 
expansion/contraction [161 1171 1181 119] . by investigating the Hamiltonian of such a system through canonical 
reduction. Several interesting exact solutions to the iV-body equilibrium problem |20l 1211 1221 123] in (1+1) 
dimensions have also been obtained. 

In this paper, we follow a similar approach in (2+1) gravity to obtain canonical equations of motion to analyze 
a two-body system. The analysis of the iV-body problem in (2+1) dimensions also has an interesting history, 
beginning with construction of a spinning point-particle solution [4] and then a consideration of the quantum 
scattering problem [24]. Further developments came upon realizing that the problem could be analyzed from a 
topological perspective [35], and an implicit solution for the metric and the motion of N interacting particles was 
obtained [26]. Based on a mapping from multivalued Minkowskian coordinates to single- valued ones, it becomes 
explicit for two particles with any speed and for any number of particles with small speed. It is possible to show 
that the collision of point particles in 2+1 AdS spacetime can result in the formation of a black hole [27] . 

The connection between these approaches and more traditional canonical methods as employed in (1+1) [H] 
and (3+1) dimensions [2H] nas n °t been explicated. In this paper, we address this issue. We begin by finding the 
total action corresponding to the system, which consists of the Einstein-Hilbert action for the field as well as a 
term corresponding to coupling gravity to matter. A variational approach will then lead to coordinate conditions 
and constraints, which will in turn produce the total Hamiltonian as an expansion in powers of the gravitational 
coupling k. The Hamiltonian will be explicitly calculated to second order in k for the general two-body case and 
to every order for the quasi-static limit. Our results agree with previous work in (2+1) dimensions which has 
explored, through geodesies, the equations of motion in the quasi-static approximation 25 . 



2 Action 

A first step is to derive the action for an iV-body system. This action will consist of two parts: the first, Ie, 
is the Einstein-Hilbert part due to the geometry of spacetime, and the second, 7m, is due to matter-gravity 
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coupling. Our derivation will follow the general ideas from the original ADM article [29] . In this formalism, the 
metric is defined as 

ds 2 = -Nodt 2 + 9H {dx 1 + N^tj (dx 3 + N 3 dt) (1) 

where No = (— g 00 ) 3 and Ni — got are the lapse function and shift covector. We also define the extrinsic 
curvature K%j = (2A r o) _1 (N^j + Nj\i — 3»j,o) and the canonical momentum conjugate to the metric: tt zj — 
— ^/g(K lJ — Kg 13 ), where the vertical bar denotes a covariant derivative with respect to the induced metric gij. 
Using the Gauss-Codazzi equations, we can perform a 2 + 1 decomposition of the Einstein-Hilbert action by 
writing 3 R as a combination of 2 R = R and extrinsic curvatures: 

Ie = -— / d z x\J~ 3 g 3 R 
2k J 

= ^ j d\N a ^g (r~K 2 + KijK ij ^j - (2^gK), + [2^g(KN' - // A,,,i| _ 

= Yk j d3xNa ^9 R ~ 2n <° + 2 ~ V99 iJ N ,j\ 4 (2) 

= 2^ / d%X {■''"■■ + NoR ° + NiRi ~ 2% -° ~ 2 l* tlN i - nN * + N j],ij . 

^hj d3x {" •" : + NoR ° + NiRi } - 

where R° = -^(gR + ty 2 — iTijir 13 ) and R l = 2ttM, and where in the last step we discarded total divergences. We 
now calculate the Lagrangian corresponding to matter-gravity coupling, which will give us the second half of the 
action. We consider the Lagrangian which minimally couples N particles to gravity: 



Cm = ]T / ds a l -m*J ( SM^-^- )\- (3) 



ds a ds a 

We can rewrite this Lagrangian in a more convenient form by introducing the canonical momentum p a(J = a ^Jf 
conjugate to z all as well as Lagrange multipliers X a which ensure PaPa/i = — m 2 : 

Cm = j ds a fzw^ 2 - - ^'ais^iPa^Pavg 1 "" (x) + ra 2 )J <5 (3) (x - Z a (So)) 

= [P^O - \ X a{Pa^Pavg^ {x) + ml) \ S {2) (r a (x°) ) (4) 



= [P aii a - N y/paiPajg ij +ml + N t p aj g 1J ) 5 (2) (r a (x )) 



where \ a — ^'a^v (because of the 5- function integration) and r a = x — z a , and where in the last step we used 
([T]) as well as the solution to the constraint PaPa^ = — wi 2 ,: 



Pao = Nip a jg 13 - No^/paiPajg* 1 ' +m%. (5) 
Combining this with our previous result for the Einstein-Hilbert action, we write the total action I = Ie + J Cm 



^ Pas i^ (2) (r a ) + ^(7r"0y,o + NoR° + NiR^ , 



(6) 



where [41129] 



R° = —={gR + 7r 2 - n lJ n lJ ) - 2k^2 \/paiPajg ij +m2,5 (2) (r a ) 

^ " (7) 

R l = 2n' l j+2Kj2p a j9' j 5 (2) (r ei ) 

a 

must both vanish following variations of No and Ni. These constraints (0 have a very natural physical in- 
terpretation. The latter of these is a momentum-balance constraint, indicating that in this (2+l)-dimensional 
self-gravitating system the particle momenta must be cancelled by the momentum of the gravitational field. The 
former is an energy-balance constraint, which basically indicates that the energy of the gravitational field (a 
combination of "potential energy" expressed as the curvature of the spacelike slice and "kinetic energy" consist- 
ing of the square of the gravitational momentum) must equal the relativistic energy of the particles. 



2 



We now consider rewriting this action by choosing convenient coordinate conditions. We define hij — gtj — <5ij, 
where <5y is the Kronecker delta, and use the orthogonal decomposition 



h 



•j — 1 £^&idj ) h + hi j + hj t i, 



where 1/ A is the inverse of the flat-space Laplacian. Note that h T and hi can be found to be [29] 



hi 



2A 



(8) 



(9) 



With these new definitions in hand, we can write 

(1 + h T )7Y"dt log(l + h T ) - ( h 



d 3 'xn 13 d t gij = j d 

= Jd 3 x -\og(l + h T )d t (-K U (l + h T ))-2w* 3 j d i (hi-^dih 



(10) 



= / d 3 xl -Alog(l + h T )<9 t 



^((i + ft T K 



where we discarded total derivatives in arriving at the last expression. Hence, the action ([6]), with constraints 
R = and i? 1 = from (Q), can be rewritten as 



1=1 d 3 xy^ paiil8(x — z a ) — ^-Alog(l + h T )do 
I z — ' Zk 



-((l+ftv*) 



-- 7 r*jd (hi-—d i h T ). (11) 



This action suggests that we should adopt the coordinate conditions 



I([l + ^« 



x = hi - ^9ih T , 



(12) 



such that, defining the Hamiltonian density H = -^A log(l + /i T ) and using x^ v — the action can be rewritten 



1 = J d 3 xl^ Pal il8{yi-7, a )-H^ 



(13) 



By imposing proper boundary conditions, we rewrite our coordinate conditions in a more convenient form: 



TV — 



_ 1 

9ij,i — ~^9jj, 



(14) 



where the second equation implies that hi — j^dih T — 0, which in turn means that hij = 8ijh T , and thus, 

g w = «i i (l + 0- ( 15 ) 

In the coordinate system defined by these coordinate conditions and defining <f> = log(l + h T ), the constraint 
equations R° = and R' — lead to 



A<^ = -2k ^(m 2 + ple-^h^-Za) - eW' 

a 

3j{'K %i e' t ') = -K)^Pai8(x - z a ), 



(16) 
(17) 



where the Hamiltonian is given by H — — ^- J d 2 xA(p and where p 2 = p 2 is the square of the norm of the vector 
p. The solution satisfying the condition n" = is 



K 

'2t7 



22 D ijk {p a )dk log i 



(18) 



where Djjfc(p a ) = PaiSjk + PajSik — PakSij- Thus, Equation (|16|) can be rewritten as 

,V2 



A<f> 



2k 2~2 (ml + pie 5(x-Zo) 

a 

j(^) Ajfc(Pa)Aji(pb)9fc log rg^j log r b j e" 



(19) 



We can then determine the Hamiltonian for a system of particles, which can be written down in the general 
form: 

fl = EH+fr* W ) 1/2 + ^ / d 2 J^]>> Ijfc (Pa)A j; (Pb)^ log r a ft log r b \e-* . (20) 

a J I a b ) 

where cj> is obtained by solving (|19[) . 

As will be exemplified in the following section, this Hamiltonian can be solved perturbatively in k. Writing 
4> — ^2, i^fi , it is clear that the first term of on the right-hand side of (|19[) only depends on <f>^ , . . . , while 
the second depends only on <j)^\ . . . , (j)^~ 2 \ Thus, we inductively solve for each term in the power series for 4> 
by calculating Atj)^' from the known <jP^ , . . . , </>' 1-1 ' , and solving it for (f>^ before moving onto the (i + 1) term. 
Provided that we can solve these A(/> differential equations, this method allows us to solve for the Hamiltonian 
up to any order in k. 



3 Two-Particle System 



To deal with the particles' divergent self-energies, which tend to be common in this sort of calculation, we 
introduce a renormalization scheme similar to that used in quantum field theory, consisting of introducing an 
energy scale to our measurement of distances. We first introduce the density of particle i to be 

Mr) = 9 Ar V 2 ( 21 ) 
iirryr + a) 

where a is an intrinsic length scale depending on the particle's energy, and where we will be taking the limit as 
a goes to zero. Notice that for any value of a, j v ndV = 1 and n(r) goes to a 8 function as a becomes small. 
Moreover, the equation Ag; = rij(r) has the solution (up to a term with vanishing Laplacian) 



91 2, 

as opposed to simply g oc log(r) were we to use the S function. Then, the self-energy of a particle will be 
log(l) = 0. a is to be thought of as the gravitational length scale of the interaction, which we assume to be on 
the order of the Planck length. Since a is dimensionful, we can take the limit of it being a very small quantity, 
while at the same time choosing units such that a — 1. Then, for a non-zero r, log(l + r/o) = log(l + r) ~ log(r). 

We will consider a system of two particles of masses mi and m.2, in the centre of inertia frame where 
pi = — P2 = p, using the notation E a = J m| + p% with, again, = p a 2 - We also use shorthand notation 
n = \x — Zi\ and r = \zi - z 2 \. 

We first note that, in this frame, 

Dijk(pa)dkh(r a ) = {pidj +Pjd t - 5ijPkdh)(h(ri) - h(r 2 )) (23) 

a 

such that the Hamiltonian (|20[) will have the form 

H = E ( m « + P' e " 0(Za) ) + ^ fd 2 * MKri) - h(r 2 )) ■ V(A(n) - fc(ra))] e~+. (24) 

a 

We will be using a perturbative approach to third order in k, and therefore calculate: 

, ,(1) . 2 ,(2) 3 ,(3) 



(m 2 a 



- < J " - \*' l '* m ) (25) 



Note that in the limit k — > 0, we wish to retrieve H = yj m\ + p%, and therefore require <f> to be of order at 



4 



least ft. The contribution to the Hamiltonian will be separately calculated for each order in k. We find: 
A0 (1) = -2 [Bim(n) + E 2 n 2 {r 2 )] 



A<f> 



(2) 



-L^^niin) + -^ (1) (z 2 )n 2 (r 2 ) 



p 



■V(ft(n) - ft(ra)) • V(A(n) - ft(ra)) 



(3) 



7T 

£1 
i_ 



Ex 



+ 1 h(r)m(n) + + 1 A(r)n 2 (r 2 ) 



e 2 



,(1)/ n2 . P i(l)f \2 
( Zl ) + 7^2^ ( Zl ) 



4?r 



A [A(n) - ft(n,)] s 



-f\ 2 



2tt 2 



fc» 1 + 



£ 2 



AE 2 

2 



ni(n) + (l«2) + f, [V(/i(n) - h(r 2 )) ■ V(ft(n) - /»(ra))] 



1 77i2t,/ \2 . P tti 2 i./ ^2 

E 2 h(r) + 



4E 2 tv 2 



E 2 h(r) 



ni(n) 



+ (1 o 2) 



^ [V(/i(n) - fc(r 2 )) ■ V(fc(n) - h(ra))] (BiA(n) + S 2 ft(r 2 )) 



where we used the solutions 



[Eih(n) + E 2 h(r 2 ) 



^(2) = z£_ 

2tT 2 



-i- I ) Mr 1 )/;(r) + ( + l j /;(;■, 



47T 2 



(ft(n) - />(r 2 )f 



(26) 



(27) 



The Hamiltonian H = / d 2 xA(j) can now be found as an expansion in 



i5i + E 2 



2tt 
2 



£2 
£1 



2 ] P_ 

4tt 2 



£2 



8tt 2 ^ 4tt 2 2 8£ 2 7r 



2^2 ^2 



(28) 



£ !og(r)' 
+ (lo2) 

+ |^(£ 1 +i?2)log(r) 2 . 

where we used equations (|39l) and (|40f) and took the limit h(r) — > log(r). In the case of equal masses, Ei = E 2 , 
our Hamiltonian is greatly simplified: 



2 up 



2 2 

K p 



H = 2E+^ log(r) + log(r) 2 (p 2 + 2£ 2 ) 

Simplifying further, the massless approximation gives: 

2k j, 3k 2 p 3 , 2 
H = 2p-\ p log r H 5— log r 



(29) 



(30) 



which we notice as being the second-order expansion of - 1 2 J r Vt^ (^plogr) where VK is the Lambert 14 7 function. 
Note that the Hamiltonian in 1+1-dim dilaton gravity was also expressed in terms of the W function |18] . 



4 Exact Quasi-Static Approximation 

We can also use this formalism to describe exactly the Hamiltonian dynamics of a system in the quasi-static 
limit. In the following, we will therefore retain all orders in the coupling constant k and only the lowest (linear) 
order in p 2 . A. Bellini et. al. treated the geodesic equations in this approximation 25 . For simplicity we 
shall consider only the equal-mass case. To simplify the notation, we will, in contrast to the previous section, 
return to a 8 mass distribution and a log(r) potential. We will explicitly keep the divergent terms logrn and 
logr 22 , where ru = \xi — Xi\, to display their cancellation at every order in k. From ([20} the Hamiltonian in this 
approximation is given by 



jj —2m + — c - mlogri2 f / g-^ mlo s''n _|_ e "i i°g >"22 
2m V 

fcp 2 f 2 

+ I d x (V(logri - logr 2 ) • V(logn - logr 2 )) 



m (log r! +log j- 2 ) 



(31) 
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where cj) a — — — logr a is a static solution to (|19[) and where mi = m-2 = m and p 1 = — p 2 = P- The integral 
of the 3rd term on the right hand side of (|31l) is 



m(log T-i+log r 2 ) 

m(log 7-i+log r 2 ) 



/" d 2 :r [V(log n — log r 2 ) ■ V(log n — log r 2 )] e 

= J d 2 a; [V(logri + logr2) • V(logri + logr2) — 4(Vlogri • Vlogr2)] e 

_ ~27T r -^m(log rn+log i- 12 ) , g * m(log n 2 +log r 2 2 ) _ 2] 

Km I i 
— 4 / d 2 x(Vlogri • Vlogr2)e 



The first term in the preceding expression cancels the second (divergent) term in eq. (|3ip . Using the integration 
formula 

d 2 x(V log n • V log r 2 )(Iog n + log r 2 ) n = — 7r(log r) n+1 , (33) 

we can rewrite the Hamiltonian (13111 as 



# = 2m + ^-^ / d 2 x(Vlogn ■ Vlogr2)e- m(logri+logr ' 2) 
m 7H ' 

P 2 



2m + P-e^ loST (34) 



m 



which indeed agrees with equation (|29[1 in the quasi-static limit. This Hamiltonian gives us the canonical 
equations of motion 

r = — = — (r ) - p (35) 
dp m 

dH 2k 2 , , , 

P = --5— = P (r ) * r, (36) 

which lead to the equations of motion in the second order 

Km 2r(r ■ r) — r(r) 2 



(37) 



This equation is identical with Eq.(4.7) in [25]: rj r = AGM , where r\ r = r] x + i r/ y ,G = f- and M = 2m 



Vr 



5 Conclusion 

We have presented here the general expression ()20p for the Hamiltonian describing N particles coupled to 
(2+1) dimensional gravity. Beginning with a derivation a general form for the action corresponding to iV-body 
dynamics, we employed the ADM formalism and proper coordinate conditions to find this general Hamiltonian. 
Our results are complementary to those that employ topological methods [251 1261 127| . 

Our result 11201) has the advantage that it can be studied in a wide variety of physical regimes, including 
small and large mass, small and large momenta, and small and large gravitational coupling. We explicitly 
demonstrated this in two cases. First we obtained an expansion to second-order in the gravitational coupling 
constant k for the two-particle system. Second we also considered a quasi-static, two-particle equal-mass system 
to every order in k, in which case our results agree with those previously found in the literature [251126] . 

A number of interesting problems remain for future consideration. An obvious thing to try is to incorporate 
additional couplings to electromagnetism and a cosmological constant. A study of the quantization of the 
Hamiltonian ([20p should also afford interesting insight into the nature of (2+1) quantum gravity coupled to 
matter. 
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6 Appendix: Fourier Integrals 

This appendix will contain proofs of some integration formulas used in section [3l which are done through Fourier 
transforms. We begin by noticing that 

zLJ d * kA ±e ik * = S( X ) (38) 
and therefore, the two-dimensional Fourier transform of log(r) is -gr-- We use this to show that, for n > 2, 

j d 2 xA(\ogr 1 . . . logr„) = — ^ J d 2 xA J f[ k ^ x e~ tT " h *' Zi 



ki e 



(39) 



(271 

= 0, 

meaning that the integral of the laplacian of any polynomial function in log r will vanish. Similarly, we can show 

d 2 x(V log n ■ V log T2) log T3 

— f d 2 xd 2 k 1 d 2 k 2 d 2 k 3 {k } \ k2 } e Kk 1+ k 2 +k 3> x e -i(k 1 .z 1+ k 2 .z 2+ k 3 .z 3 ) 



(2tt)3 J " *"* ,V1 " "" 3 kik 2 2 k 2 



3 



L f ^j, ,2, ,2, ( fc l ' feg) _ ... , n „-j(fe 1 -2 1 +fc 2 .« 2 +fe 3 ^a) 



d feid fc 2 d fc 3 „ g „ 5(ki + k 2 + k 3 ) e" 



7T ( log r 3 i log r 32 - log n 2 log n 3 - log r 2 i log r 23 ) . 



which implies that 

/ d 2 a;(Vlogri • Vlogr2)logri = — 7r(logr) 2 

(40) 



d 2 a;(Vlogri • Vlogr2)logr2 = — 7r(logr) 2 

These integrals will not change when we make the replacement log(r) — > log(l + r/a) for suitably small a, as 
they get little contribution for the small r part of the integral. This can be seen by noticing that replacing the 
bounds of integration by J Q E makes the integrals vanish. 
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